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1. INTRODUCTION 


Calculating the electronic energy band structure of solids 
by means of the linear-combination-of-atomic-orbitals one usually 
keeps interactions only up to nearest neighbours. Nevertheless, 
it is knomm, that for example, for the body centered cubic (bcc) 
crystal lattice the strength of the next-nearest-aeighbour (nnn) 
interactions is not small in conparison with the nearest-neigh- 
bours ones [1]. In [2] it has been shown that inclusion of nnn 
interactions has greatly improved the accuracy of the calculated 
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(in tight-binding method) valence band structure and density of 
states for diamond and zincblende crystals. Kiwi M. et al. [3] 
have investigated the influence of nnn interactions and overlap 
between orbitals centered on neighbouring atoms m the band 
Structure of the Bethe and diamond lattice. They showed, that 
inclusion of nnn interactions changed considerably the electronic 
structure of these solids. 

The purpose of the present paper is the calculation of the 
density of states and electrical conductivity of disordered al- 
loys with nnn interactions taken into account. The remainder of 
this paper is arranged as follows. In Sec. 2 we present a model 
alloy Hamiltonian and investigate the influence of nnn interac- 
tions on the density of states of a pure bee crystal. Next, in 
Sec. 3, we introduce the nnn bcc lattice Green functions and give 
the recurrence equation for them. In Sec. 4 we present the for- 
mulas for dc and ac electrical conductivity with nnn interaction 
included. Numerical results for a density of states of a pure 
nnn bec crystal and for a binary alloy crystal we present in 
Sec. 2 and 2, respectively, and numerical results for a de elec- 
trical conductivity in Sec. 4. 


2. THE ALLOY HAMILTONIAN 


Let us consider a substitutionally disordered binary alloy 
A,B,., with atoms A, B randomly distributed at sites of a cubic 
perfect lattice. The alloy is described by a single s-like con- 
duction band and is assumed to be treated in the tight-binding 


approximation in Wannier's representation. The Hamiltonian reads: 
H= Qin > ea EZ» ho <al = (1a) 
n m 


(nn) (nnn) 


= ŻIn »éQ«Q4 +2, 2 Lin d<ml+ h, È Sin > nl (15) 
n ném ném 
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Here the |n?» are *annier functions localized at the n-th 
site. The diagonal elements may be regarded as random atomic 
levels which take one of two possible values Si or €x depending 
on whether an atom of type A or B occupies the site n. As 
we can see from Eq. (1) we assume the hopping integrals Baan 
responsible for the propagation of electrons, to take on a non- 
-zero value only when "n" and "m" are nearest-neighbours (nn) 
and next-nearest-neighbours CEŁ In general hopping integrals 
can take six possible values 20 ni), nf), nĘ), SEI nf? 
according to the occupation of sites "n" and "m", but in the 
following we assume constant h, and ho values, independent on 
the nature of the atoms located at nearest or next-nearest-neigh- 
bour lattice sites. This means, we consider diagonal disorder, 
only. We assume the lattice of the alloy to be body centered cu- 
bic. The Hamiltonian of the ideal crystal composed, say of A-type 
atoms reads: 


(an ) (nnn ) 
H, = & Zln><nl +h, 7 Zlny<n| +h, Z$Sln><nml (2) 
n ném (ném 


This Hamiltonian in "k", Bloch-like representation defined by 


res? ics 
\k> = N 2e 13» G) 
j 


is diagonal, and the eigenvalue €(k) is given by 


gc uu 
€(k) = €, - Bh, cos — cos cos - 
A 1 2 2 2 
(4) 
ak ak 
- 4h, [ cos? = + cos? m s + coa- WK 6n, 
2 2 


So the pure A(B) component band is centered at €, (€p) and has 
a width SÉ GER = 16h, (the band extends from JULI - 8h, - 


- 6h, to €. (Eg) + 8h, - 655). 
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3. THE DENSITY OF STATES AND "nnn" LATTICE GREEN FUNCTIONS 


The lattice Green functions for Hamiltonian with only 
nearest-neighbour hopping integrals included has been widely 
used [5-12]. Because of additional interactions included in the 
Hamiltonian (1) we must introduce the "nnn" lattice Green func- 
tions. They are defined in the following way: 


Ti 
G (1m; E,a )= = V 


^ 
o  E* + cosx cos ycos Z+ a (cos^xecos yscos^z ) 


cos lx cos mx cos nx dx dy dz 


(5) 
where Et = E + io. 
The density of states of the ideal crystal may be written 
in terms of these "nnn" lattice Green functions as: 


k 
DEE) = e gi, Im 6.000; z, äer ax k= Yn, (6) 


Problem of calculations of the density of states for a pure bec 
lattice with "nnn" transfer integrals included has been investi- 
gated in papers [13-15]. Nevertheless, for the sake of complete- 
ness and for latter discussion we calculate this density of states 
for several values of the ratio of transfer integrals h, and hy 
and present in Fig. 1. For a computational convenience we used 

the half bandwidth as an energy unit. The inclusion of the "nnn" 
transfer integrals in Hamiltonian (1) changes the shape of the 
density of states in a rather clearly way. 

Firstly, the unphysical divergence in the middle of the band 
disappears. Secondly, the maximum moves in the upper energy band 
side (for h5/h4 < 1) and simultaneously decreases and broadens. 

We can.see, especially near the maximum of the curve D (Fig. 1) 
for hh, = 0.5 some similarity to the shape of the SC density 

of states (with "nn" transfer integrals in Hamiltonian included, 
only). This similarity increases with increasing ratio ho/h, aid 
in limit hy/h4 —» co one can obtain the "nnn" sc tight-binding 
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EN -0.6: ENERGY 0.6 1.0 
| 
Fig. 1. The density of states for a pure "nnn" bcc crystal (in 


tight-binding representation). The parameters ho/h are: A - 0.0, 
B — 0.1, © — 0.3, D — 0.5. Energy in units of half-bandwidth. 


density of states (for s-like states) - see, for comparison 
Di, 15]. 


For a lattice Green functions defined in Eq. 


(5) we can 
write the following equation: 


(E + Axe ET i c(i +$ ET 2 GT +8 8,5) = & 5 
(7) 


The first summation is performed over lattice vectors connecting 
nearest=neighbour atoms and the second one over the next-nearest- 
-neighbour atoms. From this equation, taking various crystal vec- 
tors 1 into consideration, we can obtain relations between dif- 
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ferent "nnn" lattice Green functions. For instance, taking 
1- (000) we have: 


Œ + 2x) G(00015,5) + G(11155,5) + ŻG(200;8,) = 1 (8) 
and for 1 = (111): 


(Œ + 2x) G (000;E,Ż) + GG (222:5,5) + 36 (220;E,% ) + 56 (002;By5) + 
+ 6(000;B,5)] + 8 x [6(111;E,5) + 6(113;E,5)]= 0 (9) 


Not all of the Green functions we get in Eqs. (8) and (9) are 
independent. In a case of the "nn" lattice Green functions it is 
known [21] that for the body-centered cubic lattice there are 
only three independent Green functions, and these are usually 
taken as G(000,E), G(200;E), G(220;E). 

The other Green functions can be obtained from them through 
relations similar to these represented in Bos, (8) and (9). In 
our case the number of independent Green functions is greater, 
but still the relations (8), (9) and those obtained from Eq. (7) 
for other = vectors may be helpful in calculations. 


4. DC AND AC ELECTRIC CONDUCTIVITY 


In the random alloy problem many investigations have been 
devoted to the calculation óf the electronic density of states, 
electric conductivity, Hall coefficient (see for example [16-19]). 
For these problems the coherent potential approximation (CPA) 
and its generalization (off diagonal disorder, cluster effects, 
short-range-order, long-range-order and so on) have been found 
useful. We use the coherent potential approximation formulated 
by Soven [20] (see also CPA method for ternary alloys [191). 

The quantity we wish to determine first for the alloy is: 


<G@)> = GE)=[E-F- 3G) (10) 
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D(E),6(E,) IN ARBITRARY UNITS 


-0.8 -04 0.0 0.4 0.8 
ENERGY, FERMI ENERGY 


Fig. 2. The density of states for the alloy described by para- 
meters: A — h/h, = 0.3, x = 0.35 B — hh, = 0.3, x = 0.5. 
Positions of an atomic levels are: A — e. EB -> G, PŁ 


B — E g >J, H (see fig. 3). 


D(E) IN ARBITRARY UNITS 


i 
-0.8 -0.4 0.0 0.4 0.8 ENERGY 
Fig. 3. The density of states for the alloy described by para- 
meters: x = 0.3 and hh, is 0.0, 0.3, 0.5 for curves A, B, C, 
respectively. Positions of atomic levels: A ebe DACH — LK; 


B — Eas Ep > GE; C — €, $5 —J ,H. 
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the configuration average of the single-electron Green function. 
Here the operator W is the second term on the right-hand side 
in Eq. (1a) and Z(E) is the self-energy operator. In a single 
site approximation (we neglect the statistical correlations be- 
tween the atomic scattering matrix and the effective waves coming 
from other sites) the CPA self-energy (E) is obtained as a so- 
lution of the following equation Dei: 


3G)2x€,Q« (1-x)ég « (€,- Z8)) 6G)(€5 - I(E)) 
(11) 
The knowledge of the self-energy enables us to calculate elec- 
tronic density of states for our model binary alloy. Results are 
shown in Figs. 2 and 3, ke have calculated the density of states 
for four values of the ratio ho/h, i.e. for 0.050015 0.5 and 
0.5. In Fig. 2 we present density of states for hj/hi = 0,5 and 
for two values of the concentration of A-type atoms: x = 0.5 
and x = 0.5. Ne can observe, that the peak's position do not de- 
pend on the concentration of various sorts of atoms, in contrast 
to the situation in the next Fig. 5. Here we have density of 
states curves for a constant concentration x = 0.3 and for dif- 
ferent values of a ratio h,/h, = 0.0, 0,5 and 0.5. As hj/h, is 
raised the peaks of the D(E) move in the upper band energy side 
and simultaneously the main peak decreases and that's one cor- 
responding to the"impurity band" increases. So, changing the 
ratio h/b] we can obtain the density of states the same as for 
"nn" Hamiltonian but for different positions of band centers and 
concentrations, 
The dc electrical conductivity tensor is given by the Kubo- 
-Greenwood formula [16]: 


2 
e . 2T eb [at =) <Tr[p" Ae sch ô Œ - E)]> , (12) 


where e, m denote the charge and the mass of electron, £2 - 
the volume of a crystal, p - the momentum operator, f(E) - the 
Fermi-Dirac distribution function and H - the Hamiltonian of the 
system. The sign <....> denotes the average Over configura- 
tions. Since this method was described elsewhere [16, 18], we 
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2 ^ 2 
2e 1 af Kë m Y = 
62a! Jaz cz d "rego 3 AP) 


where summation is carried over the Brillouin zone of the bcc 
crystal, v Y (k) is the w-th component of the electron velo- 
city and Im G(k,4*) = In[a4* - ex, - Z(4* j1^! is the spec- 
tral density of states of the averaged crystal. Here we calculate 
exactly the function Fk) contrary to commonly used approximate 
formula [16], and finally obtain 


685) = = H CET hf. (x,y,2)]- 


2 
pa e dx dydz (14) 
B, ei. Re TG)? + [In TAJ 


where: 


^ 3 
f4 (x:y,2) = sin?x cos?y cus? + cos?xsin^y cos z 4 cos?x cos^ysin^z 


2 2 


f, (x,y,2) = cosx cos y cos z (sin^x + sin?y * sin^z) 


4 (sin? 2x + sin? 2y + sin?2z) (15) 


fz (x;352) 


The frequency dependence of the real (absorptive) part of the ac 
electrical conductivity tensor is described by [16]: 


e In, Ze” Se s la-ta H e tr[v sin 5 vf SD -H)]> + 
(16) 
Following the way described, for example in [4] we obtain 


2 
= A AT — a 
Elo) = £— Re | an sip — ferm 
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Es D-O) Pto) HOLT ihan, 


E649) - ZG) Ely to- IN 
Here the function $ reads (in terms of a "nnn" bcc lattice Green 
functions): 
2 


a^h 
© (z) = za [(1 + k2)6(000;z4,k,) + 6 (200;2,,k,) - 6(220;z,,k,) - 
- G (222$24 ,k4) * 4kG (111524, Kk.) - 4k G (511524, X4) - 


2 
- k*G (40032, ,k S 
(4005 24 ,k, )]; (18) 
h h h 
3 2 GN: 2 M EZ 
s lapi: vp, eU e ce IU 


Numerical results of the dc electrical conductivity are pre- 
sented in Figs. 4 and 5. The conductivity is plotted as e func- 
tion of the Fermi energy Eg» and half of the bandwidth of the 
pure material is taken as a unit of energy. In Fig. 4 we show 6 for 
x = 0.3 and for two values of ratio ho/h, equal to 0.5 - curve A 
and 0.5 - curve B. Ve can see that conductivity depends on a va- 
lue of the next-nearest neighbour transfer integral, it is great- 
er for greater ratio Baba, This is consistent with the physical 
picture of the motion of electrons in systems described by Hamil- 
tonian (1), namely, electrons have an additional possibility for 
propagation im crystal. 

In order to give a simple and transparent picture of the in- 
fluence of "nnn" transfer integrals on dc conductivity we present 
in Fig. 5 conductivity curves for the alloy described by concen- 
tration of A atoms equal 0.5 and hj hà = 0.5 = curve A and hh, = 
= 0.0 = curve B. In a case of "nn" Hamiltonian a conductivity is 
a symmetrical curve (and density of states, too). For the "nnn" 
Hamíltonian situation is quite different. Now conductivity is a 
strongly asymmetric function of Ev and an "impurity" part of this 
curve has a smaller value in comparison vith the "host" part, al- 
though is located in energy region of increased values of a den- 


Sity of states. 
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6 (Ez) IN ARBITRARY UNITS 


-08 -04 0.0 04 0.8 


FERMI ENERGY 
Fig. 4. The dc electrical conductivity for the alloy described 
by parameters: x = 0.3 and h./h, equal to 0.3, 0.5 for curves 
A and B, respectively, 


DIE) IN ARBITRARY UNITS 


-0.8 -0.4 0.0 0.4 0.8 ENERGY 


Fig. 5. The dc electrical conductivity (curve A) and density of 
states (curve B) of the alloy described by parameters: x = 0.5, 
h,/b, = 0.5. Curve C - dc conductivity for "nn" Hamiltonian = 


fOr comparison. 
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In conclusion, we can say, that the conductivity curves as 
functions of Fermi energy level for the "nn" and "nnn" Hamilto- 
nians may be quite different depending on the ratio ho/h4. This 
means, that when we want investigate the behaviour of the de 
conductivity for a specific model Hamiltonian of a solid ihe pa- 
rameters of ihis Hamiltonian must be chosen very carefully. 


REFERENCES 


W. M'OowrpiN XN. S. „U k a dae DĄ and „JK pengo S2: t] Bhysoes Soc. 
Japan, 1974, 37, 1278. 

2. C hapati D,uJ.pand Cfo hfe n M. L.: phya., stat. solą by 
1976, 68, 405. 

3.4K 2 wi «Mr, „Rami r e2 Ra ands Tr i a a A.ż Phys. 
Rev. B, 1978, 17, 3063. 

4. Wy S o Kd ak Xvid... Ta r an k-o.. 5. and Tła ma n= 


ko R.: J, Phys. C, 1980, 13, 6659. 
D UMOSMICEC! TG ES tee Sie 152 1912,25. 165, 
6. Glasse r M. D.: J. Math. Phys., 1972, 13, 1145. 
Ne JXofyECRERNG. NSt SJENRES. AE ten D592 
8. JacozYacE (GONG: J.EPhys. 0, 1970, 4, 1510, 
9. Liao reiotwaw=Te=mand="Hec"r"i"g-u*c-h-l+M.:=Je=llath.=PhySs., 


1972, 13, 1243. 

10.Guttmann-4A.J. and Joyce G. 8.: J. Phys. A, 
EE 

MM. M or ita rege, It eme Arel Eau, 
1971, 12, 986. i 

12. ori i 8 T., Horiguchi and Chen C. C.: 
J. Phys. Soc. Japan, 1971, 38, 981, 

3. Swendsen R.H. and Callen H.: Phys. Rev. B, 
1972, 6, 2860, 

14, Lo ly P. D.: Phys. Rev. B, 1973, 8, 4405. 

15.Loly P. D. and Buchheit M.: Phys. Rev, B, 1972, 
5. 1986, 

lO. ETT Ee R. Jog K num h afue l1 m^ eA. Mand 
Lea th P. L.: Rev, Mod. Phys., 1974, 46, 465. 

TL. H o shel nol eke eeu PhYS.-Cy 1977, tO, 2147; 


Influence of the Next Neorest-neighboure-interactions ... 501 


1a. Ve Li c k $ B.: Pays. Rev., 1369, 184, 614. 

19. wysjo kad as KM WKP anen PTT a t Jj: J-Phys:C; 
1976, 9, 4271. 

20. So ven P.: Phys. Rev., 1967, 156, 809. 

21. H ori t a t.a*dnePnyeswi, 1975 mE 478. 


STRESZCZENIE 


W pracy obliczono stałoprądowe przewodnictwo elektryczne 
stopu opisanego hamiltonianem ciasnego niązania zawierajacya 
różne od zera całki przeskoku dla najbliższych i dalszych sąsia* 
dów. Wyniki obliczeń numerycznych wskazują, że dołączenie do 
opisu całek przeskoku dla dalszych sąsiadów całkowicie zmienia 
zależność przewodnictwa od energii Ferniego. 
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